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We investigate the entanglement in a two-qubit Heisenberg XYZ system with different Dzyaloshinskii- 
Moriya(DM) interaction and inhomogeneous magnetic field. It is found that the control parameters (D x , B x 
and b x ) are remarkably different with the common control parameters (D Z ,B Z and b z ) in the entanglement and 
the critical temperature, and these x-component parameters can increase the entanglement and the critical tem- 
perature more efficiently. Furthermore, we show the properties of these x-component parameters for the control 
of entanglement. In the ground state, increasing D x (spin-orbit coupling parameter) can decrease the critical 
value b xc and increase the entanglement in the revival region, and adjusting some parameters (increasing b x 
and J, decreasing B x and A) can decrease the critical value D xc to enlarge the revival region. In the thermal 
state, increasing D x can increase the revival region and the entanglement in the revival region (for T or b x ), and 
enhance the critical value B xc to make the region of high entanglement larger. Also, the entanglement and the 
revival region will increase with the decrease of B x (uniform magnetic field). In addition, small b x (nonuniform 
magnetic field) has some similar properties to D x , and with the increase of b x the entanglement also has a 
revival phenomenon, so that the entanglement can exist at higher temperature for larger b x . 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 

Entanglement is one of the most fascinating features of 
quantum mechanics and it provides a fundamental resource 
in quantum information processing (j]]. As a simple system, 
Heisenberg model is an ideal candidate for the generation and 
the manipulation of entangled states. This model has been 
used to simulate many physical systems, such as nuclear spins 
1 2], quantum dots 01, superconductor J3l and optical lattices 
1 5[], and the Heisenberg interaction alone can be used for quan- 
tum computation by suitable coding |0] . In recent years, the 
Heisenberg model, including Ising model 0], XY model Hi, 
XXX model U, XXZ model El and XYZ model JH 
have been intensively studied, especially in thermal entangle- 
ment and quantum phase transition. Very recently, F. Kheiran- 
dish et dl. and Z. N. Gurkan et al. lfl2ll discussed the Heisen- 
berg XYZ model with DM interaction (arising from spin-orbit 
coupling) and magnetic field. However, almost all the direc- 
tions of spin-orbit coupling and the external magnetic field are 
fixed in the z-axis in the above papers. The external magnetic 
field along other directions is discussed only in some spe- 
cial Heisenberg models, and the spin-orbit coupling parameter 
along other directions has never been taken into account. This 
motivates us to think about what different phenomena will ap- 
pear if we choose the parameters along different directions in 
the generalized Heisenberg XYZ model. To explore this, here 
we choose x-axis as the direction of spin-orbit coupling and 
the external magnetic field in the Heisenberg XYZ model. 

In this paper, we discuss the differences between the 
Heisenberg XYZ models with parameters in different direc- 
tions, and then analyze the influences of these parameters on 



the ground-state entanglement and thermal entanglement. We 
find that x-component DM interaction and inhomogeneous 
magnetic field are more efficient control parameters, and more 
entanglement and higher critical temperature can be obtained 
by adjusting the value of these parameters. In order to pro- 
vide a detailed analytical and numerical analysis, here we take 
concurrence as a measure of entanglement [13]. The con- 
currence C ranges from to 1, C = and C = 1 indi- 
cate the vanishing entanglement and the maximal entangle- 
ment respectively. For a mixed state p, the concurrence of 

the state is C(p) = max{2A max — Ym=i ^i)0}> where A^s 
are the positive square roots of the eigenvalues of the matrix 
R = p(a y (g)(T y )p*(a v ®a v ), and the asterisk denotes the 
complex conjugate. If the state of a system is a pure state, 
i.e. p = |*)(*|, |*) = a|00) + 6|01) + c|10) + the 
concurrence will be reduced to C(|^)) = 2\ad — bc\. 

Our paper is organized as follows. In Sec. II, we com- 
pare the DM coupling parameters along different directions. 
In Sec. Ill, we compare the external magnetic fields (includ- 
ing uniform and nonuniform magnetic field) along different 
directions. In Sec. IV, we analyze the entanglement of gener- 
alized Heisenberg XYZ model with x-component DM inter- 
action and inhomogeneous magnetic field. Finally, in Sec. V 
a discussion concludes the paper. 



H. THE COMPARISON BETWEEN THE DM COUPLING 
PARAMETERS ALONG DIFFERENT DIRECTIONS 



A. DM coupling parameter along z-axis 
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The Hamiltonian H for a two-qubit anisotropic Heisen- 
berg XYZ model with z-component DM coupling and external 



2 



magnetic field is 



H 



JxO x x al + J y o\ol + J z a\a z 2 + D z (afa^ - a\a 2 ) 



-(B z + b z )<j\ + (B z -b z )a z 2 , 



(1) 



where J,(i = x,y,z) are the real coupling coefficients, D z 
is the z-component DM coupling parameter, B z (uniform 
external magnetic field) and b z (nonuniform external mag- 
netic field) are the z-component magnetic field parameters, 
a l (i = x, y, z) are the Pauli matrices. The coupling constants 
Ji > corresponds to the antiferromagnetic case, and Ji < 
corresponds to the ferromagnetic case. This model can be re- 
duced to some special Heisenberg models by changing Jj. Pa- 
rameters Ji, D z , B z and b z are dimensionless. 

Using the similar process to lfl2ll . we can get the eigenstates 
of If: 



|0i, 2 > =sin0i, 2 |OO) +cos0i, 2 |ll) 



(2a) 



1^3.4) - sin0 3 ,4|Ol) + X COS #3,4 1 10), (2b) 

with corresponding eigenvalues: 

E 1 , 2 = J z ±wi, (3a) 



where wi 



E 3 A = ~Jz ± W 2 , 

V4i? 2 2 + (J x ~ Jy) 2 , 



(3b) 



"'2 



v/462 + ADl + (Jx + J y ) 2 , and X = 



?X, 2 = arctan( J* Jy B ), 9 3A 



J x + J v -2iD z 

^± Wi ^2b;>> "3.4 - arctan^ ± W2 ^ 2bz )• 
The system state at thermal equilibrium (thermal state) is 

p(T) = cxp{ JSt\ w here Z = Tr[exp( 7 =^)] is the partition 
function of the system, H is the system Hamiltonian, T is 
the temperature and Kb is the Boltzmann costant which we 
take equal to 1 for simplicity. Then we can get the analytical 
expressions of Z, p(T), and R, but we do not list them 
because of the complexity. After straightforward calculations, 
the positive square roots of the eigenvalues of R can be 
expressed as: 



ZW2 



^co S h 2 (^)-4&2 S inh^) 



2/^2. 



T 



T 



,w 2 



±smh(^) x / (J x + Jy) 2 +4D 2 



(4a) 



e T 



^3,4 



Zwi 



^cosh 2 (H)-4Bfsmh^(^) 



2/^1 



T 



T 



(4b) 



±sinh(^-)(J a: - J y ] 



Thus, the concurrence of this system can be written as 11311 : 

C = max{|Ai - A 3 | - A 2 - A 4 , 0}, (5) 

Fig. 1(a) demonstrates the concurrence versus temperature 
T and D z , where J x = 1, J y = 0.5, J z = 0.2, and 




(b) 




FIG. 1: (Color online) (a) Thermal concurrence versus T and D z 
without the external magnetic field, (b) Thermal concurrence versus 
T and D x without the external magnetic field. Here J x = 1, J y = 
0.5, and J z = 0.2. 



B z = b z = 0. It shows that the concurrence will decrease 
with increasing temperature T and increase with increasing 
D z for a certain temperature. The critical temperature T c is 
dependent on D z , increasing D z can increase the critical tem- 
perature above which the entanglement vanishes. These re- 
sults are in accord with those in Ref. 11211 . 



B. DM coupling parameter along x-axis 

The Hamiltonian H for a two-qubit anisotropic Heisen- 
berg XYZ model with x-component DM coupling and exter- 
nal magnetic field is 

H = J x a\a 2 + J y o\a\ + J z a z a 2 + D x (afa 2 - o\a\) 



-(B x + b x )af + (B x -b x )a%, 



(6) 



similarly, where Jj(i = x,y,z) are the real coupling coeffi- 
cients, D x is the x-component DM coupling parameter, B x 
(uniform external magnetic field) and b x (nonuniform exter- 
nal magnetic field) are the x-component magnetic field pa- 
rameters, and a % {i = x, y, z) are the Pauli matrices. The cou- 
pling constants Jj > corresponds to the antiferromagnetic 
case, and Ji < corresponds to the ferromagnetic case. This 
model can be reduced to some special Heisenberg models by 
changing Jj. Here, all the parameters are dimensionless, and 
we introduce the mean coupling parameter J and the partial 
anisotropic parameter A in the YZ-plane, where J = Jy + Jz 
and A = 

In the standard basis {|00), |01), 1 10) , |11)}, the Hamilto- 
nian H can be rewritten as: 



H 



J z 
Gi 

Jx J 11 



Jx 



G 2 

— Jz 

+ 

G a 



G 3 

J x ~\~ Jy 
Jy J z 

Go 



J x Jy 

Gi 

Jz 



(7) 



where Gi, 2 = iD x + B x ±b x , G 3 , 4 = —iD x + B x ± b x . By 
calculating, we can obtain the eigenstates of H 



\1> 



1 



1,2 



^(sin<y9i j2 |00) + cos(/3i j2 |01) 
V 2 

+ cosv?i i2 |10) + sin t^ 1;2 |H)), 



(8a) 



3 



2-x/b 2 -\-D 2 

^3,4 = arctan( — y g —, ). In the above standard basis, the 

— Jy—J z ±W 2 

density matrix p (T) has the following form: 



P(T) 



U x Q\ Q* 2 U 2 

Qi V x V 2 Q 2 

Qi V 2 Vx Qx 

U 2 Q* 2 Q\ Ux 



(10) 



FIG. 2: (Color online) (a) Thermal concurrence versus DM coupling 
parameter D z (red dotted line), D x (blue solid line), where T = 
3. (b) Thermal concurrence versus temperature T for D z = 3 (red 
dotted line) and D x = 3 (blue solid line). Here J x — l,J y — 
0.5, J z = 0.2, and B z , x = b, x = 0. 



where 

Ux,2 



2Z 



(e- 



sin (px + e 



sm ip 2 



J x u>2 „ Jx-\-u>2 „ 

±e T sin tp 3 ± e T sin ^4), 



1^3,4) = (sin (^3, 4 1 00) +x' c osV3,4|01) 

-X cos</? 3i4 |10) - sin</3 3i4 |ll)), (8b) 
with corresponding eigenvalues: 

E 12 = J x ±w 1 , (9a) 



E- 



3,4 



-Jx ± ^9 



(9b) 



where X = =^j0t> w 'i = V^B 2 + ( J y - J z ) 2 , w' a = 
^Abl + 4L> 2 + (J y + J z ) 2 , and 971,2 = arctan(- 



1 J-u — Jz^ZW, 



J 



Vx,2 = 7^p-{e T 1 cos 2 ipx+e r 1 cos 2 <p 2 
±e T cos (p3±e T cos fi), 



1 -fx + 
Qi, 3 = W (e 



T sin 93 1 cos tpi 



T Sill (^2 COS tp2 



±e T x sin <^3 cos v?3 ± e T x sin 934 cos 934 ) . 

Then the positive square roots of the eigenvalues of the matrix 
R can be obtained 



z' 



cosh( 



2w 2 
T 



(Co 



■sinh 2 (^)± 



UDi + {j y + j z y 



(C 9 



sinh 2 (^ 2 -) 



1662 , 
— r# sinh 



(11a) 



e t 



v 3,4 



cosh(— ) 



8B 



w 



2 sinh 2 (^)± 



T 



(Jy Jz) 2 



■ 1 2 / 2w 'l \ 

smh 2 (^i) 



165 



^ sinh 4 (^) 



T 



(lib) 



r 



where Z = 2[e^ cosh(^f ) + cosh(^l-)]. Thus, ac- 
cording to lfl3ll . the corresponding concurrence of this system 
can be expressed as: 



C = max-tlAi - A 3 | - A 2 - A 4 ,0}, 



(12) 



In Fig. 1 (b), we plot the concurrence C as a function of T and 
D x , where J x = 1, J y = 0.5, J z = 0.2, and B x — b x = 0. It 
shows the similar results to those in the last subsection, but if 
we compare the two figures in Fig. 1 in detail, we can easily 
find some differences between them, i.e., there is less disen- 
tanglement region in Fig. 1(b) than in Fig. 1(a), and increas- 
ing x-component parameter D x can make the entanglement 
increase more rapidly, for example, when T = 6 the concur- 
rence increases more rapidly in Fig. 1(b) than in Fig. 1(a). 



These phenomena show that x-component parameter D x has 
some more remarkable influences than the z-component pa- 
rameter D z on the thermal entanglement and the critical tem- 
perature. To demonstrate the differences more clearly, we plot 
Fig. 2, which shows the advantages of using parameter D x . In 
Fig. 2(a), we can easily find that for a certain temperature, D x 
has more entanglement than D z , and increasing D x can in- 
crease the entanglement more rapidly. Also, in Fig. 2(b), we 
can see that for the same D x and D z , D x has a higher critical 
temperature. In a word, increasing D x can enhance the entan- 
glement, and slow down the decrease of the entanglement (by 
increasing critical temperature) more efficiently than D z . So 
by changing the direction of DM coupling interaction, we can 
get a more efficient control parameter. 
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FIG. 3: (Color online) (a) Thermal concurrence versus uniform mag- 
netic field B z (red dotted line), B x (blue solid line), where T = 0.1. 
(b) Thermal concurrence versus temperature T for B z = 1 (red dot- 
ted line) and B x = 1 (blue solid line). Here J x = 1, J y = 0.8, J z — 
0.2, and D z , x = b z , x = 0. 




the concurrence decreases with increasing T for B z = 1 (Fig. 
3(b)) and b z = 1.5 (Fig. 4(b)). 

B. External magnetic field along x-axis 

In this case, we use the model in Eq. (6). According Eq. 
(12), we analyze similarly the variation of the entanglement 
and the roles of the external magnetic field along x-axis, by 
fixing other parameters. 

In Fig. 3(a), by comparing with \B Z \, we can see that, 
though \B X \ has some similar properties to \B Z \, \B X \ has a 
larger critical value and a more remarkable revival. In Fig. 
4(a), we find that the entanglement decreases more slowly 
with the increase of nonuniform magnetic field \b x \ than \b z \, 
and \b x \ has more entanglement when \b z \ = \b x \. Also, In 
Fig. 3(b), we find that B x has more entanglement than B z for 
a certain temperature when B z = B x . Fig. 4(b) shows that, 
for the same b x and b z , b x has a higher critical temperature 
and more entanglement (for a certain temperature). 

It is evident that for some fixed parameters, using x- 
component magnetic field can increase the entanglement more 
efficiently than z-component magnetic field. So we can 
change the direction of the external magnetic field to obtain 
a more efficient control parameter. 



FIG. 4: (Color online) (a) Thermal concurrence versus nonuniform 
magnetic field b z (red dotted line), b x (blue solid line), where T = 
0.2. (b) Thermal concurrence versus temperature T for b z = 1.5 (red 
dotted line) and b x = 1.5 (blue solid line). Here J x — 0.2, J y — 
0.4, J z = 1, and D z . x = B z , x = 0. 



III. THE COMPARISON BETWEEN THE EXTERNAL 
MAGNETIC FIELDS ALONG DIFFERENT DIRECTIONS 

A. External magnetic field along z-axis 

Here, we use the same model expressed in Eq. (1). Ac- 
cording to Eq. (4), we can write the expression of concur- 
rence, and then analyze the variation of the entanglement and 
the roles of the external magnetic field along z-axis, by fixing 
other parameters. 

In Fig. 3(a), the thermal concurrence as a function of the 
uniform magnetic field is plotted, where T = 0.1. It is shown 
that, for small \B Z \, the concurrence has a high entanglement 
and does not vary with the variation of \B Z |, and with increas- 
ing | B z | , the entanglement drops to zero suddenly at the criti- 
cal value of \B Z |, and then undergoes a revival before decreas- 
ing to zero. In Fig. 4(a), the thermal concurrence is plotted 
versus the nonuniform magnetic field with T = 0.2. It is ev- 
ident that increasing \b z \ will decrease the entanglement. In 
Fig. 3(b) and Fig. 4(b), the concurrence as a function of the 
temperature T is plotted for the uniform magnetic field and 
the nonuniform magnetic field, respectively. It is shown that 



IV. HEISENBERG XYZ MODEL WITH X-COMPONENT 
DM INTERACTION AND INHOMOGENEOUS MAGNETIC 
FIELD 

According to Sees. II and III, we know that x-component 
DM interaction and inhomogeneous magnetic field are more 
efficient control parameters of entanglement. So in this sec- 
tion, we will analyze the properties of parameters of the model 
(Eq. (6)) in detail. 



A. Ground state entanglement 

When T = 0, this system is in its ground state. It is easy to 
find that the ground-state energy is equal to 

E 2 = J x -w 1 , ifj x <\(w 1 -w 2 ), (13a) 



E 4 = -J x -w 2 , ifJ x >-{w 1 -w 2 ), (13b) 

and for E 2 and E t , the ground state is \i(j 2 ) and \ip4) respec- 
tively. At the critical point J x = \{w 1 — w 2 ) (i.e. E 2 = E 4 ), 
the ground state is \G) = ^fd^) + I ipi)) ( an equal super- 
position of \ip 2 ) and \ip4)). Thus, the concurrence of ground 
state can be expressed as 
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C(T = Q) = t i 



Jy J Z 

( Jy-Jz 



2Dt(Jy-J Z ){w 2 +Jy+J Z ) 



(14) 




In Fig. 5, the ground-state concurrence is plotted as a 
function of 6 K for different D x . With increasing 6^, the 
ground-state concurrence C keeps a constant initially, and 
then varies suddenly at the critical value of b x (b xc = 
\yJ{2J x - w[) 2 - {J y + J z ) 2 - 4L>2), at which the quan- 
tum phase transition occurs. In the region of b x > b xc , C 
has a revival before decreasing to zero. By increasing D x , 
the critical value b xc will decrease, and thus the revival re- 
gion will become larger. In addition, there will be more en- 
tanglement in the revival region for larger D x . The ground- 
state concurrence as a function of D x for different parameters 
is shown in Figs. 6(a)-6(d). These figures show that, with 
the increasing of D x , the concurrence keeps a constant ini- 
tially, which depends on J, A and B x , before reaching the 
critical value D xc = \ ^{2J X - w[) 2 - (J y + J z ) 2 - 4&2, 
then varies suddenly at D xc , and then undergoes a revival to 
reach its maximum value (C = 1) finally in the region of 
Dx > D xc . The critical point D xc decreases as b x and J 
increase, and increases as B x and A increase. Thus, by ad- 
justing the parameters of the system (increasing b x and J, de- 
creasing B x and A), we can decrease the critical point D xc to 
get larger revival region with more entanglement. 

B. Thermal entanglement 

As the temperature increases, the thermal fluctuation is 
introduced into the system, thus the entanglement will be 



FIG. 6: (Color online) The ground-state concurrence versus D x for 
different parameters, (a) b x = (red dotted line), 0.8 (green Dash- 
dot line), and 1.5 (blue solid line), where J = 0.5, A = 0.8, J x = 
-1, and B x = 1. (b) A = 0.2 (blue solid line), 0.8 (green Dash-dot 
line), and 1.5 (red dotted line), where J — 0.5, J x = — 1, B x = 0.8, 
and b x = 1. (c) B x = 0.5 (blue solid line), 1 (green Dash-dot line), 
and 1.5 (red dotted line), where J = 0.5, A = 0.8, J x = — 1, and 
b x = 1. (d) J = 0.2 (red dotted line), 0.8 (green Dash-dot line), 
and 1.5 (blue solid line), where A = 0.5, J x = — 1, B x = 0.8, and 
b x = 1. 



changed due to the mix of the ground state and excited states. 
To see the variation the entanglement, we use the thermal state 
p (T) to describe the system state, and the concurrence C to 
measure the entanglement. According to Eq. (12), the con- 
currence C is plotted in Fig. 7 and Fig. 8 by fixing some 
parameters. 

Fig. 7 demonstrates the influence of some parameters on 
the entanglement and critical temperature. In Fig. 7(A) 
and 7(a), when D x is small, with increasing temperature the 
concurrence decreases to zero at T c \ (the first critical tem- 
perature), and then undergoes a revival before decreasing to 
zero again at T c i (the second critical temperature). When 
D x > D xc , with increasing temperature the concurrence de- 
creases to zero at T c i- Increasing D x can increase the revival 
region (decreasing T c \ and increasing T C 2) and enhance the 
entanglement in the revival region, so the entanglement can 
exist at higher temperature for larger D x . In Fig. 7(B) and 
7(b), we see that decreasing B x can enhance the entanglement 
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T 



FIG. 7: (Color online) (A) Thermal concurrence versus D x and T, 
where B x = 3, b x = 1.5. (a) Thermal concurrence versus T for 
D x = (red dotted line), 1 (green Dash-dot line), D xc — 1.575 
(blue dashed line), 2 (black solid line). (B) Thermal concurrence 
versus B x and T, where b x = 1.5, D x = 1. (b) Thermal concur- 
rence versus T for B x = 5 (red dotted line), 3 (green Dash-dot line), 
B xc ~ 2.63 (blue dashed line), 1 (black solid line). (C) Thermal 
concurrence versus b x and T, where B x = 3, D x = 1.6. (c) Ther- 
mal concurrence versus T for b x = (red dotted line), 1 (green 
Dash-dot line), b xc ~ 1.47 (blue dashed line), 2 (black solid line). 
Here J x = 0.8, J y = 0.5, and J z = 0.2. 



and increase the revival region (decreasing T c x), so the entan- 
glement has its maximum value when B x = and T = 0. 
Fig. 7(C) and 7(c) show that small b x has some similar prop- 
erties to D x . In addition, Fig. 7(C) shows that the concurrence 
has a revival phenomenon as b x increases, so the entanglement 
can also exist at higher temperature for larger b x . 

On the whole, Fig. 7 can be divided into two regions: (i) 
The region with revival. In this region, D x < D xc , b x < b xc , 
or B x > B xc , where D xc , b xc and B xc are the critical val- 
ues of the parameters. When T < T c \, A max = A 3 and thus 
d = max{A 3 - X\ -X' 2 - A 4 , 0}. When T = T cl , X' 3 = X[ 
and thus d = 0. When T > T cX , A max = X[ and thus C' = 
max{X 1 — A 3 — A 2 — A 4 , 0}. With the increase of temperature, 
the entanglement decreases to zero at T C 2 ultimately, (ii) The 
region without revival. In this region, D x > D xc , b x > b xc , 
or B x < B xc . For arbitrary temperature, A max = A 1; so the 
concurrence is C = maxjAj^ — A 3 — A 2 — A 4 , 0}. Also, the 
entanglement vanishes finally at T C 2 with increasing tempera- 
ture. 

Fig. 8 shows the cross influence of the DM parameter and 
the magnetic field on entanglement. In Fig. 8(A) and 8(a), 
when D x is fixed, with increasing B x the entanglement is a 
constant initially, and then dropped to zero suddenly at B xc . 
For B x > B xc , a revival phenomenon occurs, and the en- 
tanglement becomes another constant. In addition, as D x in- 
creases, the critical value B xc increases, and the entanglement 




FIG. 8: (Color online) (A) Thermal concurrence versus B x and D x , 
where b x = 1.5, T = 0.5. (a) Thermal concurrence versus B x for 
D x = (red dotted line), 2 (green Dash-dot line), 5 (blue solid 
line). (B) Thermal concurrence versus b x and D x , where B x = 
3, T = 0.5. (b) Thermal concurrence versus b x for D x — (red 
dotted line), 1 (green Dash-dot line), 3 (blue solid line). Here J x = 
0.8, J v = 0.5, and J z = 0.2. 



enhances in the region of B x < B xc . In Fig. 8(B) and 8(b), 
when D x is small, with increasing b x the entanglement is also 
a constant initially, then dropped to zero suddenly at b xc , and 
then undergoes a revival. Furthermore, increasing D x can in- 
crease the revival region (decreasing b xc ), and enhance the 
entanglement in the revival region. When D x is large enough, 
with increasing b x the entanglement decreases from its maxi- 
mum value to zero. 



V. DISCUSSION 

We have investigated the entanglement in the generalized 
two-qubit Heisenberg XYZ system with different DM interac- 
tion and inhomogeneous magnetic field. By comparing with 
the common z-component parameters (D Z ,B Z and b z ), we 
find that the x-component DM interaction and inhomogeneous 
magnetic field (D X ,B X and b x ) are more efficient control pa- 
rameters for the increase of entanglement and critical temper- 
ature. Furthermore, we analyze the properties of x-component 
parameters for the control of entanglement. In the ground 
state, increasing D x can decrease the critical value b. xc to in- 
crease to revival region, and increase the entanglement in the 
revival region. The critical value D xc can also be decreased 
by increasing b x and J, or decreasing B x and A. In the ther- 
mal state, for T or b x , increasing D x can increase the revival 
region and the entanglement in the revival region, for B x , in- 
creasing D x can increase the critical value B xc to largen the 
region of high entanglement. In addition, decreasing B x can 
also increase the entanglement and the revival region. Small 
b x possesses some similar properties to D x , and with the in- 
crease of b x , the entanglement also has a revival phenomenon 
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so that the entanglement can exist at a higher temperature for 
larger b x . Our results imply that more efficient control param- 
eters can be gotten by changing the direction of parameters, 
and more entanglement and higher critical temperature can be 
obtained by adjusting the values of these parameters. Thus, by 
using the state with more entanglement as the quantum chan- 
nel, the ultimate fidelity will be increased in quantum commu- 
nication. In addition, in practice the need for low temperature 
enhances the complexity and difficulty of constructing a quan- 
tum computer, however by increasing the critical temperature 
the entanglement can exist at a higher temperature, this will 
reduce the technical difficulties for the realization of a quan- 
tum computer. 
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